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Введение 

При расчете технологических показателей разработки нефтяных 
месторождений часто используют различные приближенные методики, 
основанные, как правило, на решении задач одномерной фильтрации 
в системе галерей [1, 2]. Следует отметить, что в настоящее время развитие 
численных методов задач фильтрации и состояние вычислительной 
техники позволяет решать сложные плоские и пространственные задачи 
с учетом максимального количества различных факторов и сил, 
действующих на процесс фильтрации нефти в пласте [3-8, 14-18]. Однако 
подробный учет всех факторов, влияющих на процесс нефтедобычи 


затруднен из-за отсутствия достоверных способов их описания. 


Модель равновесной изотермической фильтрации многофазной 


жидкости в деформируемом коллекторе 


Процесс напорной изотермической фильтрации в упругом пласте 
описывается системой нелинейных дифференциальных уравнений 
относительно распределения насыщенности 5: 
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где 5 - насыщенность порового объема соответствующей фазой, причем 


при наличии только двух фаз должно выполняться соотношение: 


5н +5в =1; | (ЗУ. == 


В - коэффициент сжимаемости, 1/МПа; 
р - давление, МПа; 


т - пористость: 


К - абсолютная проницакмоть, м2; 
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Е(5) - относительная фазовая проницаемость; 
д - вязкость, Пас; 


рк (5) - функция капиллярного давления: 


26° соѕ0 


рк (5) = (5), (4) 


т 
где: с" - поверхностное натяжение на границе фаз, Па/м’; 
9 - угол смачивания на границе раздела жидкости и породы; 
Ј(5) - капиллярная функция Леверетта. 
Индексы у параметров: 

н - нефть; 

в – вода; 

с – скелет коллектора. 

Относительные фазовые проницаемости представляют собой 
однозначные безразмерные функции насыщенности. Типичная форма этих 
кривых представлена на рис.1. Значения относительных фазовых 
проницаемостей в зависимости от насыщенностей определяются 
экспериментально и обычно аппроксимируются степенными 
зависимостями. 


Общий вид кривых относительных проницаемостей [13] 
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Рис.1 


В [8] отмечается, что лучшее совпадение фактических и расчетных 


данных достигается при использовании квадратичных зависимостей: 


Ев = (5в Св): (1-а:(1—8в)) - (5) 
| | Ен = (1-58 -си)-(1-Ь-5в)° 


где св - насыщенность связанной воды; 


ср - остаточная нефтенасыщенность; 
а, Б – аппроксимационные коэффициенты, зависящие от свойств 
фильтрующихся жидкостей и скелета коллектора. 


Для решения задачи фильтрации многокомпонентной жидкости в 


пласте необходимо найти поле давления. Уравнение распределения 


ИА 


давления может быть получено из уравнений (1, 2) [8]: 
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где: Н - нефтенасыщенная мощность, м; 
и с - функция гидропроводности: 
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где у - вектор скорости фильтрации, тогда изменение насыщенности ^^ 
вытесняющей жидкости может быть определено следующим образом: 
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Численная реализация математической модели 
Полученная система уравнений может быть решена только численным 
методом. Уравнения в частных производных заменяются их конечно- 
разностными аналогами. 
Для случая прямоугольной декартовой системы координат покроем 
исследуемую область равномерной по х и у сеткой ячеек (плоский случай): 
хл =(1+12)А х, Ах>0, і=0/,...,№; (М +ПАх=Х и 


О = , 
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где: Ах. Ду Е размеры шагов сетки, 
* № М - число ячеек в направлении х, у соответственно (точка 
с координатами (і, /) совпадает с центром ячейки). 

Здесь используется «шахматная» сетка, т.е. координаты сеточных 
функций разнесены в пространстве, как показано на рис.2 [12]. Это дает 
возможность наглядно интерпретировать каждую ячейку как элемент 
объема среды, который характеризуется рассчитываемыми в его центре 


давлением Р; уэ а также насыщенностями компонентов. Знание же 


нормальной составляющей вектора скорости на границе ячейки позволяет 
непосредственно вычислять поток вещества через эту границу. 2 
При использовании указанного шаблона уравнение (6) может быть 


заменено явным дискретным аналогом следующего вида [9]: 
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Формула (11) проста для вычислений, т.к. использует значения, 


полученные на предыдущем шаге по времени. Однако в вычислительной 


практике явные формулы подобного типа практически не используются 
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всилу существенных ограничений на шаг по времени [9], например, 


в случае квадратной сетки (Ах =Ду): 


д < Ах? /4. (12) 


Сеточный шаблон для плоского случая 
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Для устранения указанного ограничения используют неявные схемы 


аппроксимации, устойчивые при любом значении шага по времени АЎ: 
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Точность второго порядка очевидна из самой конструкции схемы. 
Однако, реализация неявной схемы (13) представляет значительные 
трудности, из-за того, что на верхнем временном слое эта схема связывает 
значения искомой функции в пяти соседних узлах на двумерном и в семи 
узлах на трехмерном шаблоне. Указанные трудности не возникают, если 
пользоваться схемой переменных направлений или схемой расщепления. 

При реализации схемы переменных направлений система уравнений 


имеет следующий вид [9]: 
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отрицательный - добывающая), м/с. 

А При реализации схемы расщепления система уравнений имеет вид: 
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Преимуществом приведенных схем является то, что каждое 

У из уравнений (14, 15, 16, 18) можно реализовать с помощью трехточечных 
| прогонок по соответствующему направлению [9, 10, 11]. 

ове Скорость перетекания компонента через границы ячейки определяется 


явной разностной аппроксимацией формулы (8): + ; 
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г. і. 1+1 
Г для остальных трех границ ячейки скорости определяются аналогично. 
В соответствии с шаблоном (рис.2) в положительном направлении осей 


координат скорости компонентов через границы ячеек положительны. 


‚. Насыщенность компонента определяется из формулы (9): 
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Для решения уравнения (6) необходимо задание начальных и 
граничных условий. 

На границе залежи может быть задано условие непротекания 
(рлотные слои), условие свободного вытекания, условие внешнего подпора 


пластовыми водами. В первом случае граничное условие задается в виде: 


+ Ру+ = Рм» (21) 
Ч ре А 
т.е. давление за границей области р, равно давлению в граничной 


ячейке, что по уравнению (19) дает нулевую скорость перетекания. 
В случае свободного вытекания: | 
Рин = 2рм — руд» ° (22) 
а в случае внешнего подпора пластовыми водами за границами расчетной 
области задается их давление. 

В качестве начальных данных необходимо задать распределение 
давления и насыщенности компонентов, их вязкость, расположение и 
режимы работы скважин (дебит или забойное давление и цикличность), а 
также характеристики самого пласта (проницаемость и пористость). 

- В приведенной модели используется модель послойно и зонально 
неоднородного по проницаемости пласта (рис. 3). Считается, что 
проницаемость пласта при переходе от одной ячейке к другой изменяется 
случайным образом. В соответствии с методикой ТатНИПИнефти [1] 
неоднородность пласта может быть учтена с помощью методов и 
соотношений теории вероятности на основании средних значений 
коэффициентов их Параметры 


характеристик пласта и вариации. 


распределения рассчитываются на основании результатов фактического 


распределения параметров. : 
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При решении уравнения (6) с использованием схемы переменных 
направлений (14, 15) с граничными условиями (21) было выявлено 
появление нефизичных результатов, при расположении скважин в 
граничных ячейках расчетной области. Это связано с особенностями учета 


функции источника №, ; в выбранной расчетной схеме. 


‚ = ~ Модель пласта [1] 


Скважины: 1-добывающие, 2-нагнетательные 


Рис. 3 
Для устранения указанного недостатка — невозможности задания 
скважин на границе области — предлагается следующая модификация 


метода переменных направлений: 
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Использованием модифицированной схемы переменных направлений 
(23, 24) с граничными условиями (21) позволяет избежать появления 
нефизичных результатов при расположении скважин в граничных ячейках 


расчетной области. 
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На рис. 4 приведены результаты тестового расчета по предложенной 


модели. В качестве исходных данных были приняты пластовое 


давление р =15МПа , дебиты нагнетательных и добывающих скважин 


9=107* м? /с, проницаемость пласта А = 0.24.10? м2. 


Результаты тестового расчета 


а) поле давлений р, атм; знаком плюс помечены нагнетательные 
скважины, ромбами — добывающие скважины; 
б) распределение водонасыщенности для времени {= 6 мес.; 
в) распределение водонасыщенности для времени {= 9 мес.; 
г) распределение водонасыщенности для времени {= 12 мес. 


Рис. 4. 
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